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基于 CA-Markov 模型 的 多 时 间 跨 度 
土地 利用 变化 模拟 
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(1 新疆 农 业 大 学 草 业 与 环境 科学 学 院 ,新疆 ， 乌 鲁 木 齐 830052; 


2 ”新 疆土 壤 与 植物 生态 过 程 实验 室 ,新 疆 


乌鲁木齐 ”830052 ) 


摘 要 : 土地 利用 / 覆 被 变化 (Land Use and Cover Change, LUCC) 模拟 是 LUCC 研究 的 主要 内 容 和 


重要 手段 。 


时 间 间 隔 是 模拟 过 程 中 的 重要 参数 ,对 模拟 结果 精度 有 何 影响 ,有 待 深入 研究 。 以 新 疆 


玛 纳 斯 河流 域 典 型 绿洲 区 四 道 河 子 镇 为 例 ,基于 还 感 影像 提取 1975 .1985 、1995 .2000 .2005 、2010 
年 和 2015 年 的 土地 利用 数据 ,分 别 以 20 a.15a.10a 和 5 a 为 时 间 间 隔 构 建 CA-Markov 模型 ,模拟 


2015 年 土地 利用 结构 ,定量 探讨 时 间 间 隔 对 CA-Markov 模型 精度 的 影响 。 结 


表明 :1975 一 2015 


年 ,四 道 河 子 镇 LUCC 以 耕地 和 和 草地 为 主 ,期 间 耕 地 、 建 设 用 地 迅速 扩张 ,林地 、 草 地 和 未 利用 地 大 
幅 减 少 , 水 域 在 1985 一 2000 年 呈现 小 幅 增长 。 耕 地 的 增加 和 草地 及 林地 的 减少 是 研究 区 近 40 a 
LUCC 最 显著 的 特征 。 对 比 模拟 结果 与 实际 结果 ,时 间 间 隔 为 20 a、15 a,10 a,5 a hy TFOM 分 别 为 
70.35% ,69.18% ,76.32% 和 88.00% 。 基 于 2005—2010 年 转移 概率 的 模拟 结果 更 接近 于 2015 年 
实际 结果 ,适合 模拟 四 道 河 子 镇 未 来 的 土地 利用 变化 。 土 地 利用 模拟 应 依据 区 域 LUCC 特征 确定 


最 佳 的 时 间 间 隔 , 提 高 模拟 精度 。 
关 键 i: 
文章 编号 : 


土地 利用 /覆盖 变化 (Land Use and Cover 
Change,LUCC ) 是 全 球 变化 的 主要 组 成 ,对 全 球 气 候 
变化 .局 地 生态 系统 循环 以 及 人 类 生存 与 发 展 具有 
重要 影响 ,对 其 时 空 过程 的 量 测 和 动力 学 机 制 的 理 
解 已 经 成 为 科学 界 关 注 的 前 沿 内 容 "“] 。 土 地 利用 
过 程 模拟 则 成 为 LUCC 研究 的 主要 内 容 和 定量 探讨 
土地 利用 变化 过 程 驱动 机 制 以 及 评估 土地 利用 对 
生态 环境 影响 的 重要 工具 ”1 。 土 地 利用 变化 模拟 
方法 日 渐 趋 于 多 样 化 ,主要 包括 :经 验 一 统计 模型 ， 
如 马尔 科 夫 (Markov ) 模型 .逻辑 回归 (Logistic Re- 
gression) 模型 ” 等 ; 动力 学 模型 , 如 元 胞 自动 机 
(Cellular Automata ,CA) 模型 5 多 智能 体 (Multi-A- 
gent System ,MAS ) 模型 !9 等 ;集成 模型 ,如 CLUE-S 
模型 "等 ;综合 模型 ,如 土地 转换 模型 (Land Trans- 
formation Model) 等 中 。 这 些 方法 在 模拟 不 同时 空 
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尺度 的 LUCC 时 各 有 优势 。 其 中 ,元 胞 自动 机 一 马 
尔 科 夫 (CA-Markov) 模 型 结合 了 马尔 科 夫 链 长 时 间 
序列 预测 的 能 力 和 元 胞 自动 机 模拟 复杂 系统 空间 变 
化 的 优势 "1, 具有 模拟 动态 、 多 情景 LUCC 的 能 力 和 
较 好 的 模拟 精度 所 -2 ,被 广泛 应 用 于 土地 利用 变化 
模拟 与 预测 研究 。 

运用 CA-Markov 模型 模拟 土地 利用 变化 , 需 先 
从 两 期 历史 土地 利用 数据 中 提取 各 地 类 的 转换 规 
则 ,进而 预测 未 来 某 时 刻 的 土地 利用 状态 。 时 间 
间隔 是 模拟 过 程 中 的 重要 参数 之 一 "| , 现 有 研究 多 
将 预测 模拟 时 段 与 规则 提取 时 段 设置 为 相同 , 即 基 
于 过 去 和 年 的 土地 利用 变化 特征 预测 未 来 +1 年 
的 土地 利用 状态 。 也 有 研究 通过 对 历史 数据 中 土地 
利用 转移 概率 进行 内 插 来 预测 未 来 不 同时 间 尺 度 的 
土地 利用 变化 “1。 已 有 众多 学 者 在 不 同 区 域内 选 


4 : 国家 自然 科学 基金 项 目 (41561019) ;国家 自然 科学 基金 项 目 (41561103 ) 
作者 简介 : 靳 含 (1993 - ) , 女 ,陕西 咸阳 人 ,硕士 研究 生 ,主要 研究 了 
通讯 作者 : R982- ) , 男 , 新 疆 库尔勒 人 ,副教授 ,博士 ,人 研究 方向 为 土地 利用 变化 、 地 理 信息 系 统 和 遥感 应 用 
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择 不 同 的 时 间 跨 度 建立 基于 CA-Markov 的 土地 变化 
预测 模型 ,预测 时 间 间 隔 包 括 短 期 (1 ~10 a) 中、 中 
期 (10 ~20 a) 和 长 期 (20 a 以 上 )" ,模拟 结 
精度 均 能 满足 要 求 。 但 限于 数据 的 可 获取 性 ,大 部 
分 研究 中 未 对 模拟 过 程 中 时 间 跨 度 的 确定 进行 说 明 
或 探讨 ,不 同 的 时 间 间 隔 对 预测 结果 精度 有 何 影响 ， 
有 待 深入 研究 。 

新 疆 位 于 我 国 西北 干旱 区 ,区 域内 人 工 绿洲 面 
积 仅 占 9% ,但 却 集中 了 90% 以 上 的 人 口 与 95% 以 
上 的 社会 财富 2] ,是 干旱 区 的 生命 线 和 人 类 主要 聚 
居 区 。 玛 纳 斯 河流 域 位 于 天 山北 莫 中 段 ,准噶尔 盆 
地 西南 部 , 古 尔 班 通 古 特 沙漠 南 缘 , 从 南 到 北 呈 现 出 
典型 的 “山地 一 绿洲 一 荒漠 ?结构 ,土地 利用 程度 总 
体 呈 北部 沙漠 区 最 低 .中 部 绿洲 区 最 高 .南部 山区 居 
中 的 特点 1。 近 40 a 来 ,大 规模 的 水 土 开发 使 该 
流域 土地 利用 发 生 显著 变化 ,中 部 绿洲 区 扩张 速度 
较 快 ,成 为 新 疆 最 大 的 绿洲 农耕 区 和 我 国 第 四 大 灌 
溉 农业 区 。 但 在 人 工 绿洲 扩张 和 农业 发 展 的 同时 ， 


85° 25’ E 
T 


86° 30'E 


玛 纳 斯 河流 域 


45° 25'N 


44° 20' N 


43° 15'N 
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生态 环境 问题 也 日 益 突出 ,破坏 了 绿洲 生态 安全 和 
土地 的 可 持续 发 展 ”|。 因 此 ,研究 和 模拟 预测 玛 纳 
斯 河流 域 土地 利用 变化 ,为 管理 和 规划 该 流域 土地 
资源 提供 参考 。 鉴 于 此 ,本 研究 以 玛 纳 斯 河流 域 四 
道 河 子 镇 为 例 ,基于 多 期 遥感 影像 提取 土地 利用 变 
化 信息 ,以 20 a .15 a 10 a 和 5 a 为 时 间 间 隔 , 建 立 
CA-Markov 模型 ,模拟 2015 年 该 区 域 土地 利用 状 
态 , 并 与 实际 结果 对 比 , 定 量 研究 不 同时 间 间 隔 对 模 
拟 精 度 的 影响 。 


1 研究 区 概况 


四 道 河 子 镇 (44°75’ ~ 44° 83'N, 85° 67' ~ 
85°79" E ) 地 处 玛 纳 斯 河流 域 中 部 绿洲 区 , 古 尔 班 通 
古 特 沙漠 南 缘 , 东 与 老 沙 湾 镇 毗邻 , 西 与 农 8 师 121 
团 接壤 , 南 接 农 8 师 144 团 , 北 至 沙漠 边缘 ,总 面积 
为 242.17 km (图 1)。 该 区 域 属 大 陆 干旱 性 气候 ， 
夏季 炎热 干燥 ,冬季 寒冷 多 风 ,年 均 降水 134 mm ,年 
均 温 6. 8 C ,年 均 蒸发 量 1500 ~ 2 100 mm, ,年 均 风 


85° 40' 卫 85° 50’ E 


44°50'N 


44° 40'N 


85° 50’ E 


1 研究 区 位 置 


Fig.1 Location of the study area 
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速 1.7 ms”。 区 域 土地 利用 类 型 主要 以 耕地 和 
草地 为 主 ,种 植 作物 以 棉花 为 主 ,经 济 发 展 较 快 ” 。 


2 数据 与 研究 方法 


2.1 数据 来 源 与 处 理 

本 研究 采用 的 数据 包括 1975 年 的 Landsat-MSS 
影像 ,1985 年 1995 年 的 Landsat-TM 影像 ,2000 年 、 
2005 年 .2010 年 的 Landsat-ETM 影像 和 2015 年 的 
Landsat-OLI 影像 ,成 像 时 间 均 在 7 ~8 月 ,影像 云 量 
少 ,质量 较 高 。 基 于 ENVI 平台 对 影像 进行 辐射 校 
正 ,MSS 影像 采取 421 波段 假 彩 色 合 成 ,TM 和 ETM 
影像 采用 432 波段 假 彩 色 合 成 , OLI 影像 采用 543 
波段 合成 。 用 1 : 100 000 地 形 图 对 2015 年 的 OLI 
影像 进行 几何 校正 后 作为 参考 图 像 , 采 用 二 次 多 项 
式 转换 方程 对 其 它 影像 进行 纠正 ,平均 位 置 误差 控 
制 在 0.5 个 象 元 以 内 ,采用 双 线 性 内 插 法 将 图 像 重 
采样 到 30 m x30 m。 结 合 研究 区 实际 情况 和 影像 
地 物 特 征 , 将 研究 区 土地 利用 类 型 划分 为 耕地 
(C1)、 林 地 (C,)、 草 地 (C;)、 水 域 (C,)、 建 设 用 地 
(Cs ) 和 未 利用 地 (Ce)6 类 。 其 中 ,水 域 包括 天 然 陆 
地 水 域 和 水 利 设 施用 地 ,建设 用 地 包括 城乡 .工矿 、 
居民 用 地 。 通 过 人 机 交互 解 译 ,得 到 7 期 土地 利用 
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数据 (图 2) 。 
2.2 土地 利用 转移 概率 矩阵 

土地 利用 转移 矩阵 能 够 反映 区 域 在 某 一 时 段 内 
各 地 类 之 间 相互 转化 的 动态 过 程 .| ,转移 概率 为 该 
时 段 内 某 一 土地 利用 类 型 转换 为 其 他 类 型 (或 保持 
不 变 ) 的 似 然 概率 ,可 作为 马尔 科 夫 预测 中 的 状态 
转移 概率 ,计算 公式 如 下 : 


S, =P, XS, (1) 


t+1 


O<sP, <1 


Pa Ki Pi 


HY, Poe 1(i, = 1,2,.…,n) (2) 


式 中 :S,、5,41 为 t+1 时 刻 的 土地 利用 状态 ;P; 为 状 
态 转移 概率 和 矩阵 ;n 为 土地 利用 类 型 。 

将 1975 一 2015 年 划分 为 4 种 不 同 的 时 间 间 隔 
了, 分别 为 长 期 (1975 一 1995 一 2015 年 ,7, =20) .中 
期 (1985 一 2000 一 2015 年 , ?= 15; 1995 一 2005 一 
2015 年 ,7 = 10) 和 短期 (2005 一 2010 一 2015 年 ， 
T=5) ,基于 不 同时 间 间 隔 下 土地 利用 数据 提取 各 
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时 段 土地 利用 类 型 的 转移 概率 ,作为 模拟 的 参数 ,分 
别 模 拟 2015 年 土地 利用 状态 。 
2.3 土地 利用 动态 度 
土地 利用 动态 度 反 映 一 定时 段 内 各 土地 类 型 的 
变化 幅度 和 速度 ,有 助 于 了 解 区 域 LUCC 的 总 体态 
势 和 结构 变化 。 计 算 公式 二 如 下 : 
U, - 


Kt U1 100% 
— U, T ý 


(3) 


式 中 :K 为 研究 时 段 内 某 一 土地 利用 类 型 的 动态 度 ; 
U, AI U, 分 别 为 研究 期 始 、 末 某 一 种 土地 利用 类 型 
的 面积 ;7 为 研究 时 段 。 

2.4 CA-Markov 模型 

马尔 科 夫 模型 通过 系统 当前 的 状态 及 其 变化 概 
率 预测 该 系统 未 来 的 状态 ,其 核心 是 不 同时 期 土地 
利用 类 型 之 间 的 转移 概率 。 元 胞 自动 机 是 一 种 时 
间 、 空 间 状态 都 离散 ,空间 相互 作用 和 时 间 因 果 关 
系 都 为 局 部 的 网 络 动力 模型 ,通过 设 定 元 胞 .元 胞 空 
间 元 胞 状态 、 邻 域 和 转换 规则 等 来 控制 土地 利用 变 
化 模拟 过 程 2 。 

本 研究 利用 IDRISI 平 台 的 Markov 模块 和 CA- 
Markov 模块 ,基于 1975—1995 年 .1985 一 2000 年 、 
1995 一 2005 年 和 2005 一 2010 年 的 土地 利用 转移 
概率 ,分 别 以 1995 .2000 .2005 .2010 年 为 初始 年 ， 
以 20 a,15 a1I0a 和 5a 为 时 间 间 隔 模拟 2015 年 
的 土地 利用 状态 ,并 与 2015 年 实际 土地 利用 状态 
进行 对 比分 析 。CA-Markov 模型 中 设置 元 胞 大 小 
为 30 m x30 m, 滤 波 器 大 小 为 5 m x5 m, 循 环 次 
数 分 别 为 20 、15 10 .5。 

2.5 ”精度 检验 

2.5.1 解 译 精度 检验 ”人 研究 选取 Kappa 系数 检验 
遥感 影像 分 类 精度 ,计算 公式 如 下 : 
P, =P, 


P,-P, 


Kappa = (4) 
式 中 :P, 为 随机 情况 下 期 望 的 正确 模拟 比例 ;P, 为 
正确 模拟 的 比例 ; 忆 为 理想 分 类 情况 下 的 正确 模拟 
比例 。 如 果 两 个 土地 利用 变化 图 像 完 全 一 样 , 则 
Kappa =1; 当 Kappa=0. 75 时 ,表明 一 致 性 较 高 , 误 
差 较 小 ; 当 Kappa <0. 4 时 ,一 致 性 较 差 ,两 幅 图 相差 
BKR”, 

1975 年 .1985 年 和 1995 年 的 解 译 结果 分 别 采 
用 20 世纪 70 年 代 的 1 : 47 000 黑白 航 片 .80 年 代 


Fil 90 年 代 的 彩色 航 片 和 地 形 图 进行 检验 ;2000 年 、 
2005 年 .2010 年 .2015 年 解 译 结 果 经 野外 调查 和 
Google Earth 历史 高 空间 分 辨 率 影 像 检 验 。 本 文 依 
据 研究 区 土地 利用 类 型 的 空间 分 布 特征 ,选取 分 层 
随机 抽样 方法 ,在 各 期 土地 利用 分 类 图 上 抽取 
100 个 样 点 ,分别 与 上 述 参考 数据 进行 比较 ,并 将 解 
译 结果 中 正确 和 错误 的 样 点 数 代 入 公式 (4) ,得 到 
各 期 Kappa 系数 分 别 为 0. 78、0.76、0.77、0.76、 
0.790. 81 和 O. 84 ,满足 土地 利用 变化 研究 的 需 
Ra, 

2.5.2 模拟 结果 精度 检验 ”研究 时 段 内 ,区 域 土地 
利用 类 型 是 否 发 生变 化 是 LUCC 研究 关注 的 重点 ， 
也 是 评价 模拟 精度 的 重要 依据 。FOM ( Figure of 
Merit) 指 数 在 LUCC 研究 中 主要 用 于 针对 某 一 区 
域 、. 某 一 时 段 内 实际 发 生变 化 部 分 的 模拟 精度 进行 
评价 ,并 能 够 通过 对 实际 变化 部 分 和 模拟 变化 部 分 
进行 比较 ,定量 分 析 模 拟 结果 和 实际 结果 的 一 致 性 
构成 ,但 缺乏 对 实际 未 变化 部 分 模拟 精度 的 验 
证 。 为 评价 区 域 土地 利用 变化 模拟 的 整体 精度 ,本 
研究 在 FOM 基础 上 加 入 针对 实际 未 变化 区 域 的 模 
拟 结 果 ,采用 TFOM (Total Figure of Merit) 指数 作为 
评价 模拟 结果 的 标准 。 该 指数 为 研究 时 段 内 模拟 结 
果 和 实际 结果 一 致 的 部 分 与 总 体 模拟 结果 的 比值 ， 
取 值 范围 为 0 ~ 100% ,0 表示 实际 结果 与 模拟 结 
没有 重奏 部 分 ,100% 表示 实际 结果 与 模拟 结果 完全 

合 , 计 算 公式 如 下 : 


TFOM =(S, +3,)/ (5, +S, +S; + 


(5) 


式 中 :5, 指 某 区 域 在 某 时 段 内 土地 利用 类 型 未 发 
生变 化 , 且 模 拟 结 果 和 实际 结果 一 致 ;$, 指 实际 土 
地 利用 类 型 发 生变 化 , 且 模 拟 结果 和 实际 结果 一 
致 ;9; 指 实际 土地 利用 类 型 发 生变 化 ,但 模拟 结 
和 实际 结果 的 转换 地 类 不 一 致 ;5 指 实际 土地 利 
用 类 型 发 生变 化 ,但 模拟 结果 未 变化 ;5; 指 实际 土 
地 利用 类 型 未 发 生变 化 ,但 模拟 结果 发 生变 化 。 
以 时 间 间 隔 20 a(1975 一 1995 一 2015 年 ) X Ø, 
1995 年 作为 预测 起 始 年 ,土地 利用 状态 为 a ,模拟 
结果 和 实际 结果 也 为 a, 此 种 情况 (5S, ) 正确 ;1995 
作为 预测 起 始 年 ,土地 利用 状态 为 a, 模拟 结果 为 
pb ,实际 结果 为 a, 此 种 情况 (5,) 错 误 , 其 它 依 此 类 
推 ,具体 见 表 1。 


S, +8;) x100% 
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表 1 土地 利用 变化 模拟 结果 和 实际 结果 的 一 致 性 构成 
Tab.1 Agreement component of simulated results 


and observed results of land use change 


构成 初始 状态 ”实际 结果 “模拟 结果 备注 

Sy a a a 模拟 正确 
> a b b 
S; a b c 模拟 错误 
S4 a b a 
Ss a a b 

WŒ a b c 为 不 同 的 土地 利用 类 型 

3 结果 


3.1 土地 利用 变化 
3.1.1 1975—2015 年 土地 利用 动态 变化 特征 GE 
40 a 来 ,四 道 河 子 镇 主要 以 耕地 和 草地 为 主 ,各 时 期 
二 者 之 和 占 比 均 在 80% 之 上 ,2015 年 达到 92.93% 。 
该 区 域 土地 利用 /覆盖 在 各 时 期 呈现 出 不 同 的 变化 
特征 ,主要 表现 为 耕地 建设 用 地 的 增加 和 草地 水 
域 及 未 利用 地 的 减少 。 其 中 ,耕地 增加 幅度 较 大 ,从 
1975 年 的 42. 66% 增加 到 2015 年 的 76.60% ,建设 
用 地 则 从 1975 年 的 2.07% 增加 到 2015 年 的 4. 8% 。 
草地 减少 幅度 较 大 ,从 1975 年 的 37. 76% 减少 到 
2015 年 的 16.33% ,林地 从 1975 年 的 9.61% 减 少 到 
2005 年 的 1.58% ,水 域 和 未 利用 地 分 别 从 1975 年 
的 3.56% 和 4.34% 减 少 到 2015 4 0. 15% 和 0. 54% 
( 表 2)。 

四 道 河 子 镇 不 同时 期 内 的 各 地 类 动态 度 存在 明 
显 的 差异 。 耕 地 动态 度 变化 在 整个 时 期 内 呈现 “两 
头 高 .中 间 低 ”的 特征 ,1985 一 2000 年 最 低 ,2005 一 
2010 年 为 最 高 ;林地 的 动态 度 处 于 先 减 少 后 增加 的 
过 程 ,其 2005 一 2010 年 的 动态 度 为 各 土地 利用 / 覆 
被 类 型 中 最 高 ;草地 的 动态 度 在 1995 一 2005 年 最 


z: 


表 2， 四 道 河 子 镇 不 同时 期 土地 利用 结构 /km 
Tab.2 Land use structure of different periods in 


Sidaohezi Town / km’ 


年 份 Ci C, C; Ca Cs Ce 

1975 103.32 23.27 91.44 8.62 5.01 10.51 
1985 122.21 17.30 85. 69 2.14 5.93 8.90 
1995 135.23 13.66 72.03 2.33 8.14 10. 78 
2000 138.05 13.69 69. 20 2.16 8.53 10. 54 
2005 155.34 6.54 66.95 1.92 8.94 2.48 
2010 179.59 3.58 45.12 1.20 10. 68 2.00 
2015 185.50 3.82 39,55 0.37 11.63 1.30 


TE:C, — FHE , C3 -林地 ,Cs 一 草地 ,Cs — 7k dah, Cs = 建设 用 地 ,C6 = 
未 利用 地 


小 ,2005 一 2010 年 变化 最 大 ;水 域 的 动态 度 先 减少 
后 增加 ,变化 幅度 在 2005 一 2010 年 最 大 ;建设 用 地 
动态 度 变化 在 研究 期 初 和 期 末 最 高 ,1995 一 2005 年 
的 动态 度 变 化 最 小 ;未 利用 地 动态 度 的 幅度 为 前 两 
期 增加 、 后 两 期 减少 , 且 1995 一 2005 年 变化 幅度 最 
大 。 在 整个 研究 时 段 , 耕 地 水域 和 建设 用 地 的 变化 
幅度 在 研究 始末 两 端 较 大 ,土地 利用 变化 较 剧 烈 ; 草 
地 和 林地 在 2005 一 2010 年 变化 较 剧 烈 ; 未 利用 地 在 
1995 一 2005 年 变化 幅度 较 大 (图 3)。 

3.1.2 不 同时 间 间 隔 的 土地 利用 变化 转移 概率 
研究 基于 Markov 模型 得 到 不 同时 间 间 隔 下 各 地 类 
的 转移 面积 和 概率 ,能 够 定量 地 分 析 和 深入 了 解 各 
类 土地 利用 转移 面积 矩阵 .土地 利用 类 型 的 转移 方 
向 及 补充 来 源 , 了 解 研究 区 土地 利用 时 空 演 变 的 过 


程 。 


从 表 3 可 知 ,整个 时 期 耕地 的 净 增 面积 比例 最 
大 ,主要 由 草地 转 入 ,是 区 域 主要 的 土地 转移 类 型 。 
草地 时 递减 态势 ,从 1975 年 的 91.44 km FEE 
2015 的 39.55km .1975 一 1995 年 .1985 一 2000 


8 口 1975 一 1995 年 1985—20004F = 1995—20054F [2005—20104F 
6 Le 
4 bs 
2 上 
N 
e 0 A 
R -2 Ry 
-4 RY 
-6 RY 
网 
-8 要 
Su 耕地 林地 草地 水 域 建设 用 地 未 利用 地 
土地 利用 类 型 
图 3 1975—2010 年 四 道 河 子 镇 土地 利用 变化 动态 度 


Fig.3 Land use change dynamic degree of Sidaohezi Town from 1975 to 2010 


tó WN 


RI 不 同时 段 四 道 河 子 镇 土地 利用 概率 转移 矩阵 


Tab.3 Transition probability matrix of land use types in Sidaohezi Town during different periods 


时 间 间 隔 土地 利用 类 型 Ci C2 C3 Ca Cs Cs 
1975—1995 年 Ci 0. 946 7 0. 001 7 0.041 1 0.001 4 0. 006 7 0. 002 4 
C, 0. 320 7 0. 440 1 0. 156 5 0. 006 0 0.043 1 0.033 6 
C3 0. 253 7 0. 027 9 0.6717 0. 004 0 0.013 7 0. 029 0 
Ca 0.514 3 0. 050 4 0.187 3 0. 230 0 0. 000 0 0.018 1 
Cs 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 
Ce 0. 198 6 0. 020 3 0. 145 4 0. 000 0 0.013 3 0.622 4 
1985—2000 年 Ci 0. 910 8 0. 022 9 0. 054 6 0. 001 2 0. 005 7 0. 004 9 
C, 0.216 7 0. 537 0 0.119 2 0.003 8 0.050 9 0.072 5 
C3 0. 253 4 0. 017 8 0. 688 0 0. 003 5 0. 005 4 0. 031 8 
C, 0.014 7 0. 000 0 0. 128 7 0. 856 6 0. 000 0 0. 000 0 
Cs 0.0000 0.0000 0.0055 0.0000 0.9945 0.0000 
Ce 0. 163 4 0. 007 4 0. 125 6 0. 000 0 0.015 7 0. 687 9 
1995—2005 年 Ci 0.958 6 0. 000 7 0.034 7 0. 000 0 0. 004 1 0. 001 9 
C, 0. 332 5 0. 389 5 0. 266 0 0. 000 0 0.012 0 0. 000 0 
C3 0. 235 1 0. 008 6 0. 743 4 0. 001 1 0. 002 4 0. 009 4 
Ca 0.1366 0. 028 0 0.111 8 0.723 6 0. 000 0 0. 000 0 
Cs 0.0101 0. 004 1 0.017 2 0. 000 0 0. 968 6 0. 000 0 
Ce 0. 399 7 0. 038 7 0.404 4 0. 000 0 0.013 4 0. 143 8 
2005 一 2010 年 Ci 0.9507 0.005 8 0.0297 0.0029 0.0045 0. 006 4 
C 0. 429 3 0. 286 6 0.221 5 0. 000 0 0.053 8 0.008 8 
C3 0. 392 6 0.011 1 0. 542 8 0. 005 8 0. 035 7 0.012 0 
Cy 0.214 0 0.0165 0.522 6 0. 222 2 0.0165 0. 008 2 
C; 0. 095 2 0. 000 9 0. 000 0 0. 000 0 0. 903 9 0. 000 0 
Ce 0.631 8 0.010 1 0.331 1 0.000 0 0.010 1 0.016 9 


年 .1995 一 2005 年 和 2005—2010 年 草地 减少 面积 
分 别 为 19. 41 km? .16. 49 km? 5. 08 km? 和 22. 78 
km 。 减少 的 草地 主要 转移 为 耕地 ,反映 了 土地 需 
求 类 型 的 改变 , 且 各 时 段 的 转移 概率 分 别 为 
25.37% 25.34% .23. 51% .39. 26% 。 林 地 面积 呈 
下 降 趋 势 , 转 出 面积 主要 流入 耕地 ,40 a 间 的 转移 概 
率 在 21% ~43% 之 间 。 水 域 的 面积 先 减 后 增 再 减 ， 
在 1975 一 1995 年 .1995 一 2005 年 和 2005 一 2010 年 ， 
水 域 面积 变化 幅度 较 大 ,主要 转 为 耕地 且 转 移 状 态 
概率 分 别 为 51.43% 、13.66% .21. 40% 。 在 城市 化 
的 推动 下 ,建设 用 地 持续 扩张 ,主要 由 草地 转化 而 
来 。 未 利用 地 主要 转化 为 耕地 和 草地 ,在 整个 研究 
时 段 , 二 者 的 转移 概率 超过 30% ,其 中 2005 一 2010 
年 的 转移 概率 达到 96. 29% ,在 该 时 期 变化 剧烈 且 
达到 最 高 值 。 
3.2 土地 利用 变化 模拟 

分 别 将 1995 .2000 .2005 年 和 2010 年 作为 起 始 
年 ,基于 1975 一 1995 年 .1985 一 2000 年 1995—2005 
年 和 2005 一 2010 年 土地 利用 变化 数据 中 提取 的 4 
个 时 段 内 土地 利用 转移 概率 作为 转换 规则 ,模拟 


20 a,15 a、10 a 和 5 a 时 间 间 隔 条 件 下 四 道 河 子 镇 
2015 EW EHAA E Bij (R RR, Ry) ,并 与 
2015 年 实际 结果 (有 RR) 进行 比较 ,模拟 结果 见 图 4。 

基于 不 同时 间 间 隔 的 模拟 结果 与 2015 年 实际 
土地 利用 现状 之 间 存 在 差异 (图 4、 图 5)。R、R,、 
R, 模拟 结果 中 各 地 类 面积 与 实际 结果 差异 较 大 , HE 
地 差异 分 别 为 30.88 km? 36. 83 km? FI 16.79 km’, 
草地 为 18.55 km 18.68 km 和 18.31 km ,未 利用 
土地 为 7.87 km? .9.86 km 和 0.06 km ,林地 水域 
和 建设 用 地 的 变化 差异 较 小 。R, 模拟 结果 与 2015 
年 土地 利用 实际 结果 最 为 接近 ,各 地 类 面积 差异 分 
别 为 7.70 km 1. 28 km? ,8. 77 km 0. 70 km’, 
1.24 km 和 0.41 km 。 结 果 表 明 随 着 时 间 间 隔 的 缩 
小 ,土地 利用 模拟 精度 随 之 增加 。 其 中 ,各 土地 利用 
类 型 的 模拟 精度 存在 差异 ,耕地 和 草地 在 所 有 类 型 
中 模拟 精度 较 低 ,其 变化 幅度 差异 均 在 20 km 和 16 
km ,由 于 耕地 和 草地 受 人 类 活动 影响 较 多 ,驱动 因 
素 较 复杂 ,精度 随 之 降低 。 林 地 水域. 未 利用 地 的 
模拟 结果 较 稳 定 ,误差 比例 均 在 5% 以内。 建设 用 
地 的 模拟 精度 均 超 过 80% ,模拟 效果 较 好 。 
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图 4 2015 年 四 道 河 子 镇 实际 与 模拟 土地 利用 现状 图 
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Fig.4 Observed and simulated land use map of Sidaohezi Town in 2015 


国 耕 地 


2015 年 土地 利用 现状 与 模拟 结果 进行 比较 , 验 
证 不 同时 间 间 隔 的 模拟 精度 。 结 果 表 明 : Ri Ra, 
R, R, 的 TFOM 为 70. 35% 、69. 18% 、76. 32% 和 
87.70% 。CA-Markov 模型 可 以 用 来 模拟 预测 四 道 
河 子 镇 的 土地 利用 空间 分 布 。 此 外 ,Rs 模拟 结果 和 
实际 土地 利用 类 型 分 布 有 着 较 高 的 一 致 性 ,误差 较 
小 ,可 见 , 短 期 时 间 间 隔 更 适合 模拟 四 道 河 子 镇 的 土 


60% 80% 100% 
回 建 设 用 地 加 未 利用 地 


Eki 
5 四 道 河 子 镇 2015 年 土地 利用 真实 与 模拟 面积 比 


Fig.5 Observed and simulated area ratio of land use in Sidaohezi Town in 2015 


地 利用 变化 。 

4 种 时 间 间 隔 模 拟 结 果 (R;) 的 一 致 性 构成 (5;) 
在 数量 、 空 间 分 布 上 存在 较 大 差异 。 从 图 6 可 知 ,R， 
中 的 S AS, 分 别 为 65.40% 和 4.94%, S3 S, F S; 
分 别 为 2.97% .20. 89% 和 5. 80% ; R, PHJ S, FS, 
分 别 为 66. 30% 和 2. 89% ,S,,S, M S, 分 别 为 
2.90% 21.18% 和 6.74% ; R, PAS, FS, 分 别 为 
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K6 四 道 河 子 镇 土地 利用 一 致 性 和 不 一 致 性 模拟 图 


Fig.6 Map representation of agreement and disagreement components in Sidaohezi Town 


73.33% 和 2. 30% , S,.S, 和 Ss 分 别 占 1. 04% 、 
16. 11% 和 6. 54% ~ R, R, FUR, 的 构成 模式 相似 ， 
S, 所 占 比例 较 高 且 分 布 在 四 道 河 子 镇 的 东北 区 域 ， 
由 于 耕地 受 人 类 活动 .社会 经 济 发 展 影响 ,草地 未 
利用 地 多 转 为 耕地 ,在 预测 过 程 中 未 模拟 出 此 种 情 
况 ;R, P S, 和 5S, 分别 为 87.43% 和 0.27% ,S, 为 各 
模拟 结果 中 最 大 , S S, A Ss 分 别 为 0. 27% 、 
2.73% 和 9.30% ,S; 所 占 比例 较 大 ,提取 转换 规则 
发 展 过 程 快 于 预测 模拟 阶段 ,转移 概率 过 大 造成 模 
拟 结 果 较 差 ,如 草地 大 部 分 转 为 耕地 ,但 没有 考虑 该 
地 区 水 资源 是 否 满足 耕地 的 扩张 ,与 研究 区 现实 情 
况 存在 差异 。 


4 讨论 


(1) 绿洲 是 干旱 区 生态 系统 结构 的 核心 , 玛 纳 
斯 河流 域 人 工 绿洲 位 于 天 山北 坡 经 济 带 , 是 新 疆 绿 
洲 区 的 重要 组 成 部 分 , 而 绿洲 的 稳定 性 和 可 持续 发 
展 一 直 是 需要 解决 的 关键 科学 问题 之 一 ”。 
1975—2015 年 , 玛 纳 斯 河流 域 LUCC 发 生 了 显著 变 
化 , 尤 以 耕地 的 增加 和 林地 ,草地 的 减少 为 其 最 明显 


的 特征 。 随 着 区 域 经 济 社会 的 迅速 发 展 , 玛 纳 斯 河 
流域 的 绿洲 农业 耕作 区 扩张 活动 导致 土地 利用 发 生 
了 显著 变化 。 过 度 的 人 类 活动 ,使 得 大 量 自然 资 
源 被 不 合理 地 开采 和 利用 ,导致 玛 纳 斯 河流 域 生态 
环境 发 生 了 巨大 变化 。 其 中 ,耕地 面积 扩张 过 快 ,水 
资源 过 度 开采 .山区 林地 与 草地 面积 退化 及 土壤 污 
染 加 剧 等 不 当 的 人 类 活动 ,对 玛 纳 斯 河流 域 土地 的 
可 持续 发 展 造成 严重 的 影响 I。 因此 ,在 未 来 的 土 
地 利用 开发 过 程 中 ,要 适度 地 控制 耕地 的 无 序 扩张 ， 
合理 地 优化 配置 耕地 资源 ,促进 玛 纳 斯 河流 域 的 可 
持续 发 展 。 

(2) LUCC 是 典型 的 时 空 变化 过 程 ,时 间 间 隔 
是 基于 CA-Markov 的 LUCC 模拟 的 重要 参数 之 一 。 
现 有 大 部 分 研究 多 是 依据 现 有 遥感 数据 确定 模拟 时 
段 ,而 较 少 考虑 某 时 段 区 域 LUCC 特征 是 否 与 该 区 
域 未 来 变化 相 吻 合 。 本 研究 采用 长 期 .中 期 和 短期 
3 种 时 间 间 隔 来 模拟 区 域 LUCC ,结果 表明 时 间 间 隔 
为 5 a 的 模拟 结果 优 于 其 它 , 即 在 模拟 目标 年 确定 
的 前 提 下 , 较 短 的 时 间 间 隔 能 够 更 为 精准 地 反映 
LUCC 的 特征 与 趋势 ,模拟 精度 随 着 模拟 时 间 间 隔 
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的 增加 而 降低 ,SAMAT 等 的 研究 结果 也 证 实 了 
这 一 点 。 

模拟 时 间 间 隔 的 确定 应 考虑 区 域 发 展 过 程 及 其 
政策 。 近 40 a 间 玛 纳 斯 河流 域 土地 利用 发 生 了 较 
大 的 变化 ,不 同时 段 内 的 变化 幅度 有 所 差异 。 
1975—1995 年 ,林地 .草地 和 未 利用 地 变化 幅度 相 
对 较 小 ,期 间 该 流域 水 土 开发 规模 急剧 下 降 , 工 农业 
结构 调整 ,耕地 扩张 缓慢 ;建设 用 地 变化 幅度 较 大 ， 
1978 年 工业 发 展 迅 速 ,人 口 增 长 过 快 ,建设 用 地 随 
ZH MO! 。 该 时 期 土地 利用 变化 特征 和 1995 一 
2015 年 的 发 展 过 程 不 同 , 导 致 转换 类 型 不 一 致 , 模 
拟 精 度 较 低 。1985 一 2000 年 耕地 、 林 地 和 水 域 变 化 
幅度 在 整个 时 期 最 小 ,由 于 城市 化 进程 加 快 ,耕地 面 
积 扩张 速率 减缓 "| ;草地 和 未 利用 地 变化 幅度 较 平 
缓 ,该 时 期 流域 实行 毁 草 开荒 政策 5 ,草地 和 未 利 
用 地 相互 转换 ,变化 幅度 较 一 致 。 与 2000 一 2015 年 
变化 幅度 相 比 ,1985 一 2000 年 各 土地 利用 类 型 变化 
昼 度 整体 趋 于 缓慢 ,转换 概率 不 一 致 导 致 结果 偏 低 。 
1995 一 2005 年 ,林地 和 未 利用 地 变化 幅度 较 大 , 草 
地 和 建设 用 地 在 整个 时 期 变化 幅度 最 小 ,2000 年 草 
场面 积 逐 年 下 降 ,过度 放牧 导致 草场 严重 退化 ;以 裸 
地 开 殿 或 弃 耕 地 收复 为 主 ,耕地 弃 耕 后 主要 转移 为 
盐碱地 草地 ”1 。 耕 地 和 水 域 的 变化 幅度 较 平缓 ,由 
于 土地 利用 范围 继续 扩张 ,开始 有 计划 地 开发 水 资 
源 ,农田 灌溉 以 集约 、 内 涵 挖 潜 为 主 ”] 。 耕 地 草地 
在 该 时 期 与 2005 一 2015 年 面积 变化 的 差异 造成 模 
拟 精度 较 低 。2005 一 2010 年 耕地 、 林 地 、 草 地、 水 域 
和 建设 用 地 在 4 个 时 期 内 变化 幅度 最 大 , 由 于 
2000—2015 年 “三 生 用 地 ”政策 的 实施 造成 ,该 时 期 
玛 纳 斯 河流 域 的 生态 用 地 (林地 、 草 地 和 未 利用 地 ) 
和 生产 生活 用 地 (耕地 和 建设 用 地 ) 所 占 的 面积 较 
大 , 且 两 者 之 间 发 生 剧 烈 的 相互 转移 ,生态 用 地 减 
少 , 生 产生 活用 地 快速 增加 "1 ,造成 2005 一 2010 年 
和 2010 一 2015 年 土地 利用 转换 类 型 和 速率 较 一 致 、 
模拟 结果 精度 高 的 原因 。 基 于 长 期 和 中 期 的 转移 概 
率 ,模拟 精度 会 降低 ,并 与 实际 结果 差异 明显 。 不 同 
时 间 尺 度 土地 利用 变化 过 程 各 异 ,模拟 时 间 间 隔 的 


土地 类 型 模拟 结果 的 误差 相对 较 大 , ONAT- 
EVALDIVIESO 等 ”i 也 提 到 过 此 种 情况 。 

(3) LUCC 过 程 往往 受 自然 条 件 、 社 会 经 济 和 
人 类 活动 等 因素 共同 影响 ,部 分 因素 既 难 量化 又 具 
有 不 确定 性 “1 。 本 研究 关注 时 间 跨 度 对 CA-Mark- 
ov 模型 的 影响 ,但 驱动 土地 利用 变化 的 其 它 因素 也 
不 容 忽视 。 例 如 水 资源 是 干旱 区 土地 可 持续 发 展 的 
重要 制约 因素 ,有 限 的 水 资源 数量 和 分 布 状况 是 控 
制 景 观 荒漠 化 和 绿洲 化 进程 的 主导 因子 ,是 干旱 地 
区 最 为 活跃 的 自然 因素 ,对 土地 开发 利用 起 着 重要 
的 作用 “1。1950 年 起 , 玛 纳 斯 河流 域 水 利 工 程 修 
建 和 灌溉 方式 不 断 改进 ,绿洲 迅速 扩张 。 水 资源 
(ATRA ) CE 1987—1998 年 增长 速度 最 快 ,1998 
年 以 后 由 于 节 水 滴灌 技术 的 推广 ,耕地 面积 持续 扩 
张 ,水 域 面积 有 所 减少 。 随 着 玛 纳 斯 河流 域 耕地 
规模 的 不 断 扩大 ,对 农业 用 水 的 需求 量 也 随 之 增加 ， 
在 未 来 一 段 时 间 水 资源 将 处 于 严重 的 超 负 荷 运行 状 
Æo CA-Markov 模型 是 一 种 基于 转换 概率 的 预测 模 
型 ,没有 考虑 耕地 扩张 与 水 资源 总 量 之 间 的 联系 , 若 
仅 按 照 转 换 概率 进行 模拟 ,未 来 区 域内 耕地 将 不 断 
扩张 ,林地 和 草地 所 需 的 水 资源 将 无 法 保证 ,这 会 对 
生态 环境 产生 难以 估计 的 影响 ”1 ,这 是 不 符合 实 
际 的 。CA-Markov 模型 虽然 可 以 较 好 地 模拟 四 道 河 
子 镇 的 土地 利用 变化 过 程 ,但 在 建 模 过 程 中 如 缺少 
对 社会 .经济 等 因素 的 考虑 ,加 上 遥感 影像 空间 分 辩 
率 和 解 译 精度 的 影响 ,在 一 定 程度 上 制约 了 模型 模 
拟 结 果 的 精度 。 因 此 LUCC 模拟 还 要 结合 研究 区 的 
实际 情况 和 限制 因素 等 ,才能 进一步 提高 模拟 精度 ， 
得 到 更 为 接近 实际 的 结果 。 


5 结论 


研究 以 玛 纳 斯 河流 域 四 道 河 子 镇 为 例 ,基于 
1975 .1985 ,1995 .2000 .2005 .2010 .2015 年 7 期 土地 
利用 数据 ,运用 马尔 科 夫 转移 矩阵 和 土地 利用 变化 
动态 度 模型 ,对 新 疆 典 型 绿洲 区 玛 纳 斯 河流 域 四 道 
河 子 镇 近 40 a 土地 利用 变化 过 程 进行 定量 分 析 , 并 


确定 应 以 区 域 土地 利用 的 变化 幅度 .速度 和 研究 区 
的 开发 力度 ,发展 政策 等 为 基准 。 此 外 ,不同 土地 利 
用 类 型 的 模拟 精度 会 有 所 差别 ,部 分 土地 利用 类 型 
的 模拟 精度 较 高 ,如 林地 、 建 设 用 地 和 水 域 等 ,这 
类 土地 变化 幅度 相对 较 稳定 ,而 耕地 和 草地 的 变化 
很 复杂 , 受 人 类 活动 .政策 等 因素 影响 较 大 ,所 以 该 


以 不 同时 间 间 隔 构 建 CA-Markov 模型 ,模拟 2015 年 
土地 利用 结构 ,经 与 实际 数据 对 比 定量 探讨 时 间 间 
隔 对 CA-Markov 模型 精度 的 影响 。 结 论 如 下 : 

(1) 近 40 a 来 ,四 道 河 子 镇 的 耕地 面积 大 幅 增 
加 ,1975 一 2015 年 增加 了 82.18 km ,在 所 有 土地 利 
用 类 型 中 变化 最 为 显著 。 林 地 水域 和 未 利用 地 经 
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历 了 “减少 一 增加 一 减少 ”的 变化 过 程 , 净 面 积 减 
少 。 草 地 呈现 持续 减少 的 态势 ,而 建设 用 地 持续 增 
加 。 整 个 研究 时 段 ,耕地 和 草地 始终 是 研究 区 所 占 
比例 最 大 的 两 种 土地 利用 类 型 ,上 且 二 者 之 间 的 相互 
转换 较 频 繁 。 

(2) 研究 以 20 a.15 a10a 和 5 a 的 时 间 跨 度 
提取 土地 利用 转移 概率 ,分 别 构建 CA-Markov 模型 
模拟 2015 年 四 道 河 子 镇 土地 利用 变化 。 结 果 显 示 
随 着 时 间 间 隔 的 缩小 ,模拟 精度 逐渐 增加 ,基于 
2005 一 2010 年 的 土地 利用 转移 概率 模拟 的 2015 年 
土地 利用 结果 更 接近 于 实际 情况 。 

(3) 基于 CA-Markov 模型 ,通过 定量 分 析 ,探讨 
了 土地 利用 变化 模拟 过 程 中 时 间 间 隔 的 选择 (20 a、 
15 a 10 a 5 a) 对 于 模拟 精度 的 影响 ,为 土地 利用 变 
化 模拟 中 时 间 间 隔 的 选择 提供 了 一 定 的 科学 依据 。 
土地 政策 在 社会 发 展 的 各 阶段 是 有 差异 的 ,不 同 政 
策 会 导致 土地 利用 方式 发 生 不 同 的 转换 变化 。 
此 ,在 今后 的 土地 利用 预测 建 模 中 ,时间 跨度 的 选择 
还 应 考虑 该 地 区 的 历史 发 展 过 程 ,依据 不 同时 期 的 
开发 特征 ,参照 研究 区 的 现实 状况 和 限制 因素 等 , 确 
定 最 佳 的 模拟 时 间 间 隔 , 以 提高 区 域 土地 利用 变化 
模拟 的 精度 。 
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Abstract: The study of Land Use/Cover Change (LUCC) and its dynamics is crucial for human existence and 
sustainable development. LUCC simulation is the main content and effective method in LUCC research and time 
span is one of the most important parameters in the simulation process ,and thus its effect on the simulation accuracy 
needs to be further studied. In order to explore LUCC simulation accuracy at different time spans ,we took Sidaohezi 
Town of Manas River Basin in Xinjiang, China as the study area. Seven land use maps were created from remote 
sensing images (MSS images in the 1975 ,TM images in 1985 and 1995 ,ETM images in 2000 ,2005 and 2010, and 
OLI images in 2015) and the land surfaces were classified into six land cover categories; cultivated land, woodland, 
grassland , water area, construction land and unused land. Then the land use transition probability during four differ- 
ent periods ( 1975—1995 , 1985—2000 , 1995—2005 and 2005—2010) were used to establish the CA-Markov 
model to simulate land use in 2015 using long-term (20-years ) , medium-term (15-and 10-years) and short-term 
(5-years) intervals respectively. The simulation process was completed by the CA-Markov module in IDRISI soft- 
ware. Finally ,four simulated land use maps for 2015 (R,,R,,R;,R,) were compared with the real land use map in 
2015 (R) ,and the Total Figure of Merit (TFOM) indexes were used to evaluate the model’ accuracy. The results 
showed that from 1975 to 2015 , cultivated land and grassland were the dominant land use types in Sidaohezi Town. 
During the study period , cultivated land and construction land expanded rapidly , whereas forest land, grassland and 
unused land have been greatly reduced, and the water area has experienced a slight increase from 1985 to 2000. 
Throughout the study period, the increase of cultivated land and decrease of grassland and woodland were the most 
significant features of LUCC. Compared the simulated results with real land use map in 2015 , five components of a- 
greement and disagreement were detected as S| ,S,,5;,5, and S;. The TFOM indices of 20-,15-,10- and 5-years in- 
tervals were 70.35% ,69. 18% ,76.32% and 88.00% respectively. It was indicated that smaller time intervals pro- 
duced higher accuracy. The time span of 5-years was the optimum time interval to simulate the land use change of 
Sidaohezi Town from 1975 to 2015. This study provided a scientific basis for the selection of time interval in land 
use change simulation. Moreover ,the simulation accuracy depends not merely on the land use transition probability 
within a certain period of time ,but the local land policies and social development characteristics in the same time 
should also be considered. 
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